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1. Introduction 

An important goal of early vision is the computation of a representation of the 
visible surfaces in an image, in particular the determination of the orientation of 
those surfaces as defined by their local surface normals [Marr 1982, Brady 1982]. 
Many processes contribute to achieving this goal, stercopsis and structure from 
motion being currently the most studied in image understanding. Three other 
important contributing processes are shape- from-contour, shape-from-texture-- 
gradients, and shape-from-shading. Several psychophysical demonstrations show 
that shape-from-contour is significantly more powerful than shape-from-texture- 
gradicnts [Clark et. al. 1956, Gruber and Clark 1956, Braunstein and Payne, 
1969]. Similarly, Barrow and Tenenbaum [1981, Figure 1.3 ff] suggest that 
shape-from-contour is a more effective clue to shape than shape-from-shading. 

In this paper we consider the computation of shape-from-contour. Figure 1 
shows a number of shapes that are typically perceived as images of surfaces 
which are oriented out of the picture plane. The corresponding (sets of) surface 
normals (up to tilt reversal) are shown to the right. It may be supposed (see 
for example [Gregory 1973, pages 168ff]) that the slant judgements in Figure 
1 are largely determined by familiarity with regular shapes such as circles and 
squares. Figure 2 strains that hypodiesis (though Gregory proposes that we are 
/ ! *** K familiar with the shape of puddles) and suggests that die computation is based 

on more general knowledge of shapes and surfaces. The mcdiod we propose 
is based on such general knowledge, namely a preference for symmetric, or at 
least compact, surfaces. Note that the contour docs not need to be closed in 
order to be interpreted as oriented out of the image plane. Finally, Figure 
3 shows that, in general, contours arc interpreted as curved three-dimensional 
surfaces. 

We develop an extrcmum principle for determining three-dimensional surface 
orientation from a two-dimensional contour. Initially, we work out the 
extrcmum principle for the cases illustrated in Figures 1 and 2, that is, assuming 
a priori that the contour is closed and diat the interpreted surface is planar. 
Later, we discuss how to extend our approach to open contours and how to 
interpret contours as curved surfaces as shown in Figure 3. 

The extremum principle maximizes a familiar measure of the compactness or 
symmetry of an oriented surface, namely the ratio of the area to the square of 
the perimeter. It is shown that this measure is at the heart of the maximum 
likelihood approach to shape-from-contour developed by Witkin [1981] and 
Davis, Janos, and Dunn [1982]. The maximum likelihood approach has had 
some success interpreting irregularly shaped objects. The method is ineffective, 
however, when the distribution of image tangents is not random, as is the 
case, for example, when the image is a regular shape, such as an ellipse or 
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Figure 1 

Two-dimensional contours that are often interpreted as planes that are oriented with 
respect to the image plane. The commonly judged slant is shown next to each shape. 







Figure 2 

Some unfamiliar shapes that are also interpreted as planes that are oriented with respect 
to the image plane. The shape on the left is from [Gregory 1973, fig. 10.9], the others 
from [Witkin 1980, p 29 and 94]. 
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a parallelogram. Our cxtrcrhum principle interprets regular figures correctly. 
We show that the maximum likelihood method approximates the extremum 
principle for irregular figures; but that the maximum likelihood method docs 
not compute the correct slant for an ellipse. Witkin [1981, Figure 5] provides 
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Figure 3 

Some shapes that |r| Interpreted as curved three-dimensional surfaces. 
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empirical evidence that the maximum likelihood method computes a good 
approximation to the perceived tilt but underestimates the slant We prove in 
Appendix A that the maximum likelihood method consistently overestimates the 
slant of an ellipse. A more 'thorough investigation of the difference between 
the Extremum Principle and the Maximum Likelihood method is needed. 
One class of figures that are readily perceived as lying in a plane other than 
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the image plane arc skew symmetries, which are two-dimensional linear (affine) 
transformations of real symmetries. Kanade [1981, page 424] has suggested a 
method for determining the three-dimensional orientation of skew-symmetric 
figures, under the "heuristic assumption" that such figures are interpreted as 
oriented real symmetries. \tye' prove' that our extremum principle necessarily 
interprets skew ^ffifhetries aV oriented real symmetries, thus dispensing with 
the need for any^guristic assumption 'to that effect. Kanade shows that there 
is a one-parameter family of possible orientations of a skew-symmetric figure, 
forming a hyperbola in gradient space. He suggests that the minimum slant 
member of the one-parameter family ;s perceived. In the special case of a real 
symmetry, Kana|es suggestion implies that symmetric shapes are perceived as 
lying in the image plane, that is having zero slant. It is clear from the ellipse 
in Figure 1 mat'lils is not "correct. Our method interprets real symmetries 
correctly. 

First, we review the maximum likelihood method. In Section 3, we discuss 
several previous extremum principles and justify our choice of the compactness 



measure. In Section 4, we derive the mathematics necessary to extremize 
the compactness measure, and relate the extrcmum principle to the maximum 
likelihood method. In Section 5, we investigate Kanade's work on skew 
symmetry. One approach to extending the extremum principle to interpret 
curved surfaces, such as that shown in Figure 3, is sketched in Section 6. In 
the final section, we relate tliis work to the psychophysical literature on slant 
estimation and image understanding work on shape from texture. 

2. The Sampling Approach 

Witkin [1981] has treated the determination of shape-from-contour as a problem 
of signal detection. Recently, Davis, Janos, and Dunn [1982] have corrected 
some of Witkin's mathematics and proposed two efficient algorithms to compute 
the orientation of a planar surface from an image contour. Witkin's approach 
uses a geometric mode! of (orthographic) projection and a statistical model of 
(a) the distribution of surfaces in space (statistics of the universe) and (b) of 
the distribution of tangents to the image contour. We shall adopt the geometric 
model, but dispense with the statistical model in favor of an extremum principle. 

First, the geometric model. Assume that the image plane is horizontal with 

coordinates (x, y) (see Figure 4) . To obtain a plane with slant a and tilt 

r, we rotate (x, y) by r in the image plane and then rotate the image plane 

/""N by a about the new y axis. Assume that the coordinate frame in the plane 

(a, r) is chosen so that it projects into (x, y). (In Section 4, we describe this 
transformation more precisely, see also Davis, Janos, and Dunn [1982, p3].) 

Now suppose that a curve is drawn in the plane (a, r) and denote by /3 the 
angle that the tangent makes at a typical point on the curve. Let a be the 
tangent angle in the image plane at the point corresponding to /?. Then a and 
are related by: 

M\ tan . „ . 

a _ T \ — H, (2.1) 

cos a 

We now turn to the statistical model, which consists of two assumptions 
called isotropy and independence. Isotropy reasonably supposes that all surface 
orientations are equally likely to occur in nature and that tangents to surface 
curves are equally likely in all directions. More succinctly, it is assumed that the 
quantities (0, o, r) are randomly distributed and their joint probability density 
function ("density") D(0,a, r) is given by [Davis, Janos, and Dunn, 1982]: 

D(p,a,T)=-=Bino (2.2) 

it' 

We assume that the ranges of the angles are: < a < f , < r < 7r, < 
3 < n. Similarly, the density of (a, r) is given by : 
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Figure 4 

The geometry of orthographic projection. See text for details of notation. 



D[o, t) = —sin a 

It 



(2.3) 



The independence assumption requires mat the image tangents {0^:1 < i < n} 
are statistically independent. That is, it is assumed that the tangent directions 
at different points on the image curve are independent. This is only true if 
the contour is highly irregularly shaped, or if the number of samples is small. 
In any case, the assumption of independence is an inherent weakness of the 
sampling approach (see for example [Witkin 1981, p. 36]). 

Witkin [1981, p. 25 - 26] shows that the joint density function D(a, a, t) is 
given by: 



1 



sin a cos a 



D(a, a, t) = 2 

v cos 2 (ct: — r) -f- cos 2 ersin 2 (a — r) 

For the conditional density D(a\a, r) we find: 



(2.4) 
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D[a\a,T)=- 



cosa 



n cos 2 (o! — t) -j- cos 2 asm (a — r) 



(2.5) 



Denote the sample [ai,ot2,. . .a n ) by A (the sample is independent by 
assumption). It has conditional density 
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D{A\o,t)= [J D[a t \o,T) (2.6) 



By Baycs' formula wc obtain 



D[a r\A) = _m°,T)Dj*,T) 



(2.7) 



Observe that the numerator is independent of a and r. The sampling approach 
takes a random sample A and defines the most likely orientation of the plane 
{o,r) to be that which extremizes D(a, t\A). Witkin [1981] quantizes a and 
r, and describes an algorithm to find the maximizing {oj, r^). Davis, Janos, 
and Dunn [1982] develop a more efficient algorithm that first estimates a and 
r and then uses those estimates in a Newton iterative process. They provide 
evidence that their method is more accurate than Witkin's. Curiously, however, 
they state [Davis, Janos, and Dunn 1982, p 24] that "the iterative algorithm was 
not used [in the experiments they report] because the intitial estimates (whose 
computation is trivial) are very accurate and the iterative scheme often failed 
to converge to the solution". 

3. Extrcmum Principles. 

Brady and Horn [1983] survey the use of extremum principles in image 
understanding. The choice of performance index or measure to be extremized, 
and the class of functions over which the extremization takes place, are justified 
by appealing to a model of the geometry or photometry of image forming and 
constraints such as smoothness. For example, the use of extremum principles 
in surface reconstruction is based upon surface consistency theorems [Grimson 
1981, 1982, 1983, Yuille 1983] and a thin plate model of visual surfaces [Brady 
and Horn 1983, Terzopoulos 1983]. 

There are several plausible measures of a curve that might be extremized in 
order to compute shape-from-contour. First, / K?ds, where k is the curvature 
of the contour, has been investigated as a curve of least energy for interpolating 
across gaps in plane curves [Horn 1981]. Here we seek a measure of a curve 
that is extremized when the plane containing the curve is slanted and tilted 
appropriately. Contrary to what appears to be a popular belief, given an ellipse 
in the image plane, / K 2 ds is not extremized in the plane that transforms the 
ellipse into a circle. Appendix B contains a proof of this assertion. Since 
ellipses are normally perceived as slanted circles, we reject the square curvature 
as a suitable measure. 

Another possible measure is proposed by Barrow and Tenenbaum [1981, p89]. 
Assuming planarity (the torsion r is zero), it reduces to 
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ds. 



One objection to tliis measure is that it involves high-order derivatives of the 
curve. This means it is overly dependent on small scale behaviour. Consider, 
for example, a curve which is circular except for a small kink. The circular 
part of die curve will contribute a tiny proportion to the integral even when 
the plane containing the curve is rotated. The kink, on the other hand, will 
contribute an arbitrarily large proportion and so will dominate the integral no 
matter how small it is compared with the rest of the curve. This is clearly 
undesirable. For example, it suggests that the measure will be highly sensitive 
to noise in the position and orientation of the points forming the contour. 

A second objection to the measure proposed by Barrow and Tenenbaum is 
that it is minimized by, and hence has an intrinsic preference for, straight lines, 
for which dic/ds zero. This means that the measure has a bias towards planes 
that correspond to the (non-general) side-on viewing position. These planes are 
perpendicular to the image plane and have slant tt/2. 

We base our choice of measure on the following observations. 

1. Contours that are the projection of curves in planes with 
large slant are most effective for eliciting a three-dimensional 
interpretation. 

(f~*S 2. A curve is foreshorted by projection by the cosine of the slant 

angle in the tilt direction, and not at all in the orthogonal direction. 

We conclude that three-dimensional interpretations are most readily elicited 
for shapes that are highly elongated in one direction. Another way to express 
this idea is that die image contour has large aspect ratio or is radially asymmetric. 
The measure we suggest will pick out the plane orientation for which the curve 
is most compact and most radially symmetric. Specifically our measure is 

M= ,J Area) „ . (3.1) 

{Perimeter) 1 

This is a scale invariant number characterizing the curve. For all possible curves 
it is maximized by the most compact one, a circle. By compact, we mean most 
radially symmetric. This gives the measure an upper bound of 1/47T. Its lower 
bound is clearly zero and it is achieved for a straight line. It follows that our 
measure has a built-in prejudice against side-on views for which the slant is 
ff/2. 

In general, given a contour, our extrcmum principle will choose the orientation 
in which the deprojected contour maximizes M. For example an ellipse is 
interpreted as a slanted circle. The tilt angle is given by the minor axis of the 
ellipse. It is also straightforward to show that a parallelogram corresponds to a 
rotated square. Appendix C discusses die interpretation of several simple shapes. 
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In particular, an ellipse is interpreted as a slanted circle, a parallelogram as a 
slanted square, and a triangle as a slanted equilateral triangle. In Section 5 we 
extend the parallelogram result to the more general case of skewed symmetry. 

We note that the quantity M is commonly used in pattern recognition and 
industrial vision systems [Agin 1980, Pavlidis 1977, Ballard and Brown 1982] 
as a feature that measures the compactness of an object. Furthermore, we can 
show that the measure M defined in Eq. (3.1) is at the heart of the geometric 
model in the maximum likelihood approach. 

From Section 2, we see that the maximum likelihood approach maximizes 
the product of a number of terms of form 

, , COSQ! 

f{a) = =^=- . (3.2) 

cos 2 (a — r) -f- cos 2 asin 2 (a — t) 

Differentiating Eq. 2.1 with respect to the arc length sj along the image curve 
and s R along the rotated curve respectively we obtain 

Kidsi _ 1 - 3 3 > 

kr dsR f[a) 

where «:/ and /cr are the curvature at corresponding points of the image contour 
and its deprojection in the rotated plane respectively. In fact, kj = da/dsj and 
k r = dp/dsR. There is no a or r dependence in the numerator of equation 
(3.3). We can write each term k ds as (ds ds)/(p ds) where p is the radius of 
curvature. Now observe that (pds)/(dsds) is just a local computation of area 
divided by perimeter squared! Hence maximizing each f(a) in the maximum 
likelihood approach is equivalent to locally maximizing area over perimeter 
squared. In section (4) we will examine this connection more rigorously. 

Finally we note that the area, as well as the perimeter, can be obtained by 
an integral round the contour. If n is the normal to the curve then it is a 
straightforward application of Stokes' Theorem to show that 

{Area)n = - f r X dt, (3.5) 

where {Area) is a scalar quantity, and r is a vector coordinate system in the 
plane of the figure. This formula simplifies the calculations and means the 
perimeter and the area can be computed simultaneously. 

4. Extremi/ing the Measure 

We now write down the measure for a curve with arbitrary orientation and then 
extrcmizc with respect to the orientation. Let the unit normals to the image 
plane and the rotated plane be k and n respectively. The slant o of the rotated 
plane is given by the scalar product 



/*"**! 



cos a = k • n (4-1)- 

Let Tr and Tj be the contour in the rotated and image planes. A vector r in 
the image plane satisfies r ■ k = 0, and is the projection of a vector v in the 
rotated plane that satisfies v • n = 0. The projection relationship between v and 
its image r is defined by: 

r=kX (vXk) = v-(v-k)k (4.2) 

_ n X (r X k) _ (n-k)r-(n-r)k 
V ~ (n-k) - MO W> 

where X denotes vector product. Now Fa and Tj have (vector) areas Ar and 
A/ given by 



A * - I I v X dy (4.4) 

l JYr 

KI== \i, Txdt (4 - 5) 

Observe that the area vectors have the same direction as the nonnal to the plane 
containing the area. In particular, Ar is in direction n and A/ is in direction 
k. Substituting Eq. 4.3 into Eq. 4.4 and using Eq. 4.1, we find 



a _ v | kx(nXA/) 

Afl = A/ H — — 

(n-k) 

__( k ' A /). 



(4.6) 



-n 



coscr 



It follows that 



l|A*|| = E^ W 



Pf 
coscr 

(We recall that the range of a guarantees that cos a is positive. 
The perimeter lengths Pr and Pi are given by 



Pr=£ INI (4-8) 

Pi = f Vt l|dr|| (4.9) 



Substituting Eq. 4.3 into Eq. 4.8 gives 



W> 2 + |:ff «•»> 
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In general there is no simple relationship between the perimeters analogous to 
Eq. 4.6 between the areas. Nevertheless, by Eq. 4.7, 

1M = JM_ (4 . u) 

Pr Pr cos ° 

and so our extremum principle is equivalent to extrcmizing cos 2 aP R , which 
we write as 

'=l„H* 2 +wf (4 - i2) 

We extremize this with respect to the orientation n of the rotated plane, 
maintaining the constraint that n is a unit vector by a Lagrange multiplier A/4. 
This gives 



;/ f(n-k)d 



2 x ( n-dr) 2 V~* (n-k) 2 rfr 2 -(n-rir) 2 

(n-k) 2 

2\-lu_W„.JA^ ( 4 - 13 ) 



r 2 + 



(n-k)J "" (n-k)' 



+ An = 

Taking scalar products of Eq (4.13) with k and n, respectively, gives 

= (n • k)A 

(n-k)/ (n-k) 



+fr-*+K£f i *-"trJ'-* <-> 



= A 



\2\ h» . 1^2 j,2 1 /„ j r \2 



We use Eq. 4.14 to remove the integral coefficient of k in Eq. 4.13, allowing 
us to express Eq. 4.13 as a sum of the second integral and k X (k X n) times 
A. We now use Eq. 4.15 to eliminate A from this new form of Eq. 4.13. We 
recall that the unit tangent t is defined by 

- * 

ds' 

where s is arc length along the contour. It follows that 

dx = tds. 



r^ 



n 



Recalling tliat (n • k) = cos a, we find 

2J> |cos 2 a + (n-t) 2 | (t • n)dr = -k X (k X a) f Los 2 o + {n ■ t) 2 V ds 

(4.16) 
where t = dr/||dr|| is the unit tangent to the image contour. 

Let the unit vectors in the x and y directions in the image plane be i and j 
and the normal to the image plane k. By definition, cos a = k • k. The tilt r 
is defined by 



(4.17) 
(4.18) 



COST 


= 


l 


■ n 


sin a 


sinr 


— 


J 


• n 



sincr 
The tangent vector t and the normal n can be written: 



t = cos a\ -4- sin a\ (4.19) 

n = sin a cos ri + sin a sin r\ + cosak (4.20) 

where a is the tangent angle in the image. We now form the scalar products 
f~\ of Eq. 4.16 with i and j to obtain (with appropriate cancelling) 



and 



2 j> (cos 2 a -\- sin 2 crcos 2 (a — r)) * cos(q — r) cos ads 
= cos t j> (cos 2 a -4- sin 2 crcos 2 (a — r))?ds 

2 <P (cos 2 a -f- sin 2 crcos 2 (a — r))~ s cos(a — r) sin ads 
= sin r j> (cos 2 a -j- sin 2 acos 2 (a — r))^ds 



(4.21) 



(4.22) 



we can rewrite diese equations, after multiplying by factors of ±cosr and 
± sin r, as: 

2 j> (cos 2 a -f sin 2 acos 2 (a — r)) — * cos 2 (a — r)ds 

= (cos 2 <7 + sin 2 ccos 2 (a — r))?ds (4.23) 

2 i> (cos 2 a -f- sin 2 acos 2 (a — r))~ * cos(a — r) sin(a — r)<is 

= (4.24) 
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We use I-lqs. 4.23 and 4,24 in Appendix D to determine the extremizing 
orientation for a skew symmetry. We can rewrite these equations in the form: 

■-— ( — f [cos 2 o + sin 2 a cos 2 [a — t))Us) — (4.25) 



ha 

— ( f (cos 2 a -f sin 2 cxos 2 (a — r))hds ) — (4.26) 

dr ^cosia J > 

to emphasize that they correspond to extremizing with respect to a and r. 

We can implement the cxtremum method directly from Eqns. (4.25) and 
(4.26) by quantizing the tangent angle a;, slant a 3 -, and tilt Tk, replacing the 
integral by a sum. This gives good results, even though we use fixed point 
integers in our edge finder. A multi-level search speeds the algorithm by a 
factor of ten. For large slant the ratio of the greatest to least value of the 
expression is large, and the result is numerically well-conditioned. For smaller 
slants (less than about 45 degrees in the case of an ellipse) the ratio is small 
and the result poorly conditioned, so that round-off errors can be significant. 

To conclude this Section, we show that these equations arc similar, though not 
identical, to those obtained by the maximum likelihood method in the limit as 
the number of sampled tangents tends to infinity. To see this we recall from Eq. 
2.7 that this method involves extremizing D(A\a, r) with respect to a and r. 
Since the denominator is independent of a and r, this amounts to extremizing 
D(A\a, r)D[o, r). This is the same as extremizing log D(A\a, r)D(a, r). Using 
2.3, 2.5 and 2.6 we obtain: 



E = n log cos a A- log sin a — ^ log(cos 2 (aj — r) + cos 2 asin 2 (ai — r)). 

(4.27) 
where we have ignored factors of ir which will vanish on differentiation. 
Dividing E by n and taking the limit as n tends to infinity gives: 

F = log cos a j dt — j> log(cos 2 (a — t) + cos 2 tfsin 2 (a — r))dr. (4.28) 
Using the identity: 

cos 2 (a — r) -f cos 2 <7sin 2 (a — r) = cos 2 a + sin 2 <7Cos 2 (a — r) (4.29) 
gives 

F = log cos a j> dx — j log(cos 2 a + sin 2 <7cos 2 (a — r))dr. (4.30) 
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This formula is similar to Eqs 4,25 and 4.26. Thus we expect the Hxtremum 

Method to give similar results to the Sampling Method when the contour is 

sufficiently irregular. However, we can show formally that the Sampling Method 

and the F.xtrcmum Method are not equivalent. In Appendix A we show that the 

Sampling Method overestimates the slant of an ellipse. The precise discrepancy 

between die mcUiods, and its practical consequence for computing shape from 

contour is currently under investigation. 

5. Skew Symmetry 

We now consider a more general class of shapes for which the maximum 
likelihood approach is not effective. Kanade [1981, sec. 6.2] has introduced 
skew symmetries, which are two-dimensional linear (afflne) transformations of 
real symmetries. There is a bijective correspondence between skew symmetries 
and images of symmetric shapes that lie in planes oriented to the image plane. 
Kanade proposes die heuristic assumption that a skew symmetry is interpreted 
as an oriented real symmetry, and he considers the problem of computing the 
slant and tilt of the oriented plane. 

Denote the angles between the x-axis of the image and die images of the 

symmetry axis and an axis orthogonal to it (the skewed transverse axis) by a and 

P respectively. The orthogonality of the symmetry and transverse axes enable 

/""^ ; one constraint on the orientation of the plane to be derived. Kanade uses 

gradient space (p, q) [Horn, 1977, Brady, 1982] to represent surface orientations. 
He shows [Kanade 1981, p. 425] that die heuristic assumption is equivalent to 
requiring the gradient (p, q) of the oriented plane to lie on die hyperbola 

p 2 cog 2 («=i) _ g f sin 2 (a = fl = _ cos(a _ fl (5.1) 

where 

Pl==pcos l"2-J- f9Sin (~T"> ,,„ 

qx = — psral — - — l + gcosl — — J 

Kanade [1981, p. 426] further proposes that the vertices of the hyperbola, 
which correspond to die least slanted orientation, are chosen within this one- 
parameter family. This proposal is in accordance with a heuristic observation of 
Stevens [1980]. In die special case that the skew symmetry is a real symmetry, 
that is in the case that a — /? = ±7r/2, the hyperbola reduces to a pair of 
orthogonal lines [Kanade 1981, page 426] passing through the origin. In such 
cases the slant is zero. In other words, Kanade's proposal predicts that real 
symmetries are inevitably interpreted as lying in die image plane, and hence 
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Figure 5 

A. partition of a skew-symmetric object into regions for the piecewise analysis presented 
in this section. 



having zero slant. Inspection of Figure 1 shows that this is not the case. A 
(symmetric) ellipse is typically perceived as a slanted circle, particularly if the 
major and minor axes do not line up with the horizontal and vertical. 

Although Kanade's minimum slant proposal does not seem to be correct, 
there is evidence (for example [Stevens 1980]) for Kanade's assumption that 
skew symmetries arc interpreted as real symmetries. We now show that the 
assumption can in fact be deduced from our Extremum Principle. As a corollary, 
in Appendix D we determine the slant and tilt of any given skew-symmetric 
figure; only in special cases does it correspond to the minimum slant member 
of Kanade's one-parameter family. 

Instead of using Eqs. 4.25 and 4.26, we prove the result from first principles 
since this enables us to use symmetry directly. First, we partition the image 
into regions as shown in Figure 5. A typical region is shown in Figure 6a and 
we skew it by an angle a to get Figure 6b. We show that the ratio of the area 
to the perimeter squared is greatest for Figure 6a with a = 0. We show this 
is also true for the sum of the regions and take the limit as the width of the 
regions goes to zero to complete the proof. 

From Figure 6a we calculate the area A = ||A|| and the the portion Pi of 
the perimeter of the shape to be 



Ai = 2al + [b - a)l (5.3) 



f~\ 
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and 



p % = 2\JP + [b - a) 2 . (5.4) 

From Figure 6b we find 

A\ = 2a/ cos a -\- (b — a)/ cos a 
= Ai cos a 



(5.5) 



and 



P\ = Jl 2 cos 2 a -\- (b — o + / sin a) 2 -f- v// 2 cos a 2 -\-(b — a — /sin a) 2 . 

(5.6) 
We need to show that 

cos a 1 . 
< — . 5.7 

pl2 pi v ' 

i i 

This is equivalent to showing 



(2 cos a — l)(/ 2 + (6 — a) 2 ) 



f~\ < ll 2 + {b — of + 2(6 - a)/ sin a}*!/ 2 + {b - a) 2 - 2(6 — a)/ sin a V. 

(5.8) 
This condition clearly holds for 7r/2 > a > 7r/3. Assume therefore that 
< a < 7r/3. Squaring both sides of (5.8), we see that the condition is 
equivalent to 

/ 4 + (b - a) 4 + l 2 (b - a) 2 (2 + —= ™^-— ) > 0. (5.9) 

cos'' a — cos a 

On completing tlie square we see this always holds provided 

4 cos a — 1 — 3cos 2 a > (5.10) 

This is so provided \ < cos a < 1, that is for < a < |. It follows that the 
ratio of die area to the perimeter squared is maximized for each region when 
a = 0. 

We now partition the shape into n regions as shown in Figure 5. Let the i th 
block have area A\ and perimeter P\ when the skew angle is a, and denote 
the area by Ai and the perimeter by P,- when a is zero. It follows from the 
above results that 



A\=Ai cos a (5.11) 



f^ 



16 



/■"•S 



-./- 




(a) 



(b) 



Figure 6 

a. A basic region before skewing it. b. The result of skewing the region shown in (a) 
through angle a. 



rs 



and 



K . Ai 



-L2L <■ l^L 
P? Pi 2 ' 

We conclude from Eqs. (5.11) and (5.12) that 

P'i 2 > P^ cos a. 
Hence we obtain 



(5.12) 



(5.13) 



It follows that 



£P?>cos a £P t 2 

n. n 

J2 A = cos a 2 yl». 

i=l t=l 



(5.14) 
(5.15) 

v-n p/2 ^ ^n p 2' ( 5 * 16 ) 

2_,i=i ^t z^i=i -n 
and taking the limit as n tends to infinity shows that our extremum principle 
interprets a skew symmetry as a oriented real symmetry. 



2Jv=iA- . EiLi-^-t 
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6. Strategy and Sensitivity 

In Section 3, we suggested that radially asymmetric, elongated, or non-compact, 
contours most readily elicit three-dimensional interpretations. We then proposed 
a compactness measure A/P 2 , and suggested that perceived surface orientation 
corresponds to the slant and tilt that maximize the measure. We briefly discussed 
the role of cxtrcmum principles in computer vision (Brady and Horn [1983] 
present a fuller discussion), and criticized Barrow and Tenenbaum's measure on 
the grounds that it is highly sensitive to noise. In Sections 4 and 5 we analyzed 
the extremum principle and showed that it corresponded closely to human 
perception of oriented planar contours, and that it interprets skew symmetries 
as oriented real symmetries. 

In this Section we return to the discussion in Section 3, concentrating on two 
additional questions concerning our measure. First, we ask whether there are 
specific reasons why a visual system should adopt the strategy of extremizing 
our measure. We suggest that the extremum principle not only determines the 
orientation of the viewed curve, but provides an estimate of the stability of the 
interpretation. We relate this to the idea of general viewpoint. Second, we 
show that the compactness measure is relatively insensitive to noise and to the 
scale at which the image is sampled. To this end we consider two cases: (i) 
adding a sinusoid to the contour, (ii) assuming that the image can be modelled 
in terms of fractals [Mandelbrot, 1982]. 

Mewing position and stability 

Consider viewing a given planar curve from the hemisphere of all possible 
directions. Consider further the way the image changes when die viewpoint 
is shifted slightly. A smooth curve will hardly change with a slight change of 
viewpoint from most viewing directions. We call these viewing directions stable 
viewpoints. Stable viewpoints can be grouped into regions whose boundaries 
correspond to viewpoints where the contour changes rapidly for a slight change 
of viewing position. These stable regions will be quite large for images of 
smooth planar curves and smooth curved surfaces. We suggest that image 
contours are intcrpretated as curves that are viewed from stable viewpoints. 
This is the essence of die general viewpoint constraint in computer vision. We 
are able to estimate how stable a given viewpoint is. 

One way to find stable viewing positions is to define a similarity measure for 
viewpoints and then find the cxtrcma of tliis measure. Sutherland [personal 
communication] has proposed A/P 2 as a similarity measure, based on empirical 
studies of animal perception. If tliis is correct, extremizing A/P 2 corresponds 
to finding the stable viewing positions of a contour, and, as a result, die stable 
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interpretations of a contour, defined to be those which observers at general 
positions are most likely to sec. This suggests that assuming general viewing 
position corresponds to assuming interpretations around which our compactness 
measure is extrcmized and where it varies slowly. 

Wc illustrate this idea by considering an ellipse. For an ellipse the stable 
viewing position corresponds to interpreting die ellipse as an oriented circle. 
Even if the circle is slanted by as much as thirty degrees, the value of A/P 2 
hardly changes, as our algorithm demonstrates, and so the circle is a stable 
interpretation. Jf the algorithm is applied to an ellipse with large eccentricity, 
corresponding to large slant, A/P 2 changes rapidly for slight changes to the 
viewing position. So if we look at an ellipse from all possible viewing angles 
it is most likely to appear as a circle and so interpreting an ellipse as a slanted 
circle is a good strategy. 

These results are preliminary, and they will be rigorously developed in a 
further paper. 

Sensitivity of the compactness measure to noise 

The second question concerns the sensitivity of the measure to noise and to 
the scale at which an image is sampled. Since our measure involves fewer 
derivatives of the contour than measures such as that proposed by Barrow 
and Tenenbaum, it should be relatively insensitive to noise. Although our 
present algorithm and the sampling mediod are essentially equally sensitive, it 
is possible to develop a less sensitive algorithm to implement our measure. The 
sampling method inevitably involves calculating tangents, however, and hence 
is inherently sensitive to noise and scale. For example if the image contour 
is continuous, but not differentiable, the sampling method cannot (strictly) be 
applied. Similarly if die orientation of the tangent to the contour has a large 
high frequency component, the sampling approach will be very sensitive to 
noise. 

Consider the sensitivity to noise of the extremum principle. Clearly, the 
perimeter P depends on die scale at which the image is viewed, while the area 
A is much less dependent. It follows that the ratio A/P 2 varies with scale. 
If a sinusoid is added to a smooth contour, for example, the area will remain 
approximately the same but the perimeter will change significantly. This does 
not imply, however, that the extrema of A/P 2 vary with scale. For example, it 
is easy to show that if a sinusoid is added to the contour, the extremum of the 
area over die perimeter squared is effectively unchanged. 

Similarly, suppose that the image can be approximated by a fractal [Mandelbrot, 
1982]. If we measure the image at a small fractal scale constant /, the perimeter 
will be given by 
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P(l) = F/ (1 ~ d ). (6.1) 

where F is independent of /, and d > 1 is the fractal dimension. Provided / is 
small, the value of A is essentially independent of/, since the fractal tends to a 
limit contour, with finite area, as / goes to zero. Now we extrcmize A/P 2 for 
all contours which can he projected into the image contour. This is equivalent 
to extrcmizing log(A/P 2 ) which, using (6.1), we write as 

Mpa) = log A- 2 log F - 2(1 - d)logl. (6.2) 

We now extremize this over a and r (Eqs. 4.24 and 4.25), and note that the 
term involving / is independent of a and r and disappear. It follows that the 
extrema of A/P 2 are independent of the scale of viewing for an image that 
can be approximated by a fractal. 

7. Interpreting image contours as curved surfaces 

Figure 3 shows a number of contours that are interpreted as curved surfaces. 
In this section we discuss one method for extending our extremum principle 
to this general case. The key observation, as it was for Witkin [1981], is that 
our method can be applied locally, To do this, we assume diat the surface is 
locally planar. At the surface boundary, corresponding to the deprojection of 
the image contour, the binormal coincides with the surface normal. The idea 
is to compute a local estimate of the surface normal by die extremum principle 
described in the previous sections and then to use an algorithm, such as that 
developed by Terzopoulos [1983], to interpolate the surface orientation in the 
interior of the surface. The method is closely related to that proposed by Brady 
and Grimson [1981] for perceiving subjective surfaces. 

The main question concerns how to apply the extremum principle locally. 
We are currently investigating the following approach. Consider the circle of 
curvature to the space curve that is the deprojection of the contour. One way 
to define the circle of curvature is as die best fitting local circle through three 
points (or more if one assumes noisy data and makes a least squares estimate) 
on the space curve near the point in question. The circle of curvature projects 
into an ellipse. We compute die best fitting ellipse at each point on the image 
contour, and compute from it a local estimate of the surface orientation by 
finding the slant and tilt of the corresponding circle, interpreted as the circle 
of curvature. It is easy to show that so long as the surface foreshortening is 
diffcrentiable, this is a good estimator of the circle of curvature and hence of 
the local surface normal. 

It requires 5 parameters to define an arbitrary ellipse. Computing the 
best fitting ellipse in die general case is a complex nonlinear problem best 
approached using a numerical descent method, though several algorithms have 
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been published recently for computing best fitting ellipses [Book stein 1979, Agin 
1981, Nakagawa and Rosen fcld 1979, Sampson 1982]. If we assume that the 
normal to the space curve is not significantly foreshortened, we can compute 
the ellipse center and major axis from the curvature of the contour at the point 
in question. The slope of the contour at that point also defines the orientation 
of the ellipse, leaving a much simpler one-parameter problem. 

We note that perceptually the strongest local cues to surface orientation 
correspond to points of high curvature. This is consistent with our method of 
locally estimating surface orientation by fitting local ellipses.. Recall that our 
compactness measure was inspired by the observation, for example on ellipses, 
that large slant is an effective cue to surface orientation and in the case of 
an ellipse tliis produces points of high curvature. In fact, we can show that 
numerical conditioning of the estimator of slant increases monotonically with 
the slant. Conversely, for straight line portions of a contour, the curvature 
is zero, and the surface is locally planar. Hence surface orientation does not 
change along the length of the straight portion. 

One issue that remains to be studied is the interface to the surface reconstruc- 
tion algorithm. Consider the image of a triangle shown in Figure 1. There 
are, in general, three different perceived orientations of the triangle correspond- 
ing to propagating the interpretation of each of the (high curvature) corners. 
Adjacent corners give inconsistent information, and so it seems necessary to use 
a labelling approach such as that proposed by Zucker, Hummel, and Roscnfcld 
[1977]. 

7. Related work 

Gibson [1950] has argued that surface orientation is directly determined by 
certain "higher order variables" in the proximal stimulus array. What exactly 
constitutes a "higher order variable" has been the subject of extensive debate. 
Nevertheless, Gibson and his followers have proposed several such, especially 
for optical flow and texture gradients. Flock [1964, Flock et. al. 1967] have 
adopted a Gibsonian perspective on the judgement of slant in texture gradients. 
He has introduced a "higher order variable" optical slant (actually "optical 
theta" in the original) as a possible basis for monocular slant perception [Flock 
1964, Eq. 5] or, at least, as a discriminant for planarity. Flock's work implicitly 
assumes planarity [Flock 1964, p. 381], and it assumes that the tilt direction 
has been computed previously. In fact, there are many assumptions in Flock's 
paper (Freeman [1965, p. 502] counts 13). Moreover, several studies [Clark 
et. al. 1956, Gruber and Clark 1956, Freeman 1965, Braunstcin and Payne 
1969] have shown that the determination of surface orientation from texture 
gradients, for example using optical slant or as in Ikeuchi [1983], is less effective 
than its determination from bounding contour. 
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(a) 



(b) 



Figure 7 

a. A box-like figure typical of those studied by Attneave and Frost [1969]. b. A sheared 

box. 
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Attneave [1972, see also Attneave and Frost 1969] describes an approach 
to determining surface orientation that has similarities to that developed here. 
First, he argues for a Praegnanz theory of perception, which stresses economy 
principles in perception. This has meant different things to different researchers. 
The Gestalt psychologists explicitly noted the link between Praegnanz and 
minimization principles (for example the soap bubble) in mathematical physics. 
Attneave [1972, p. 285] suggests that such minima may be computed by "hill- 
climbing techniques". The extremum principles discussed in Section 3 of this 
paper and surveyed in [Brady and Horn 1983, Grimson 1981, 1982, 1983, and 
Terzopoulos 1983] can be considered to be more sophisticated formalizations of 
similar intuitions. 

In fact, Attneave adopts Hochberg's formulation of the Praegnanz theory as 
the tendency to keep differences to a minimum. In particular he considers 
three-dimensional interpretations that equalize one or more of angles, lengths 
of edges, and surface slopes in figures such as Figure 7a. The more of these 
that are in fact equalized in a particular three-dimensional interpretation of an 
image, the more likely that interpretation is to be chosen. In fact, the extremum 
principle developed in this paper will interpret Figure 7a correctly. It will also 
determine the shear in Figure 7b. 

Second, Attneave [1972, Fig. 4] considers the judged orientation of rhombus 
figures (sec Appendix D below for analysis of this case by our extremum 
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principle.) Surface orientation can be judged reliably for such shapes, more 
reliably in fact than for the box figures shown in Figure 7, so long as the 
symmetries of the rhombus do not align with the horizontal and vertical. It is 
well known that the horizontal and vertical are important for shape description, 
(see for example Brady [1983] for discussion). Note that surface slant is 
consistently underestimated. 

In order to express depth constraints, Barnard [1982] and Ikeuchi [1983] 
have proposed alternative projections to orthographic projection used in the 
development in Sections 3 and 4. The methods we have developed extend 
straightforwardly to these alternative projections. 

Barnard [1982] finds surface orientation from vanishing points of families of 
parallel edges using central or perspective projection. He considers the projec- 
tion of an angle, and, after the fashion of Hochberg and Attneave, heuristically 
assumes equi-angularity to arrive at a three-dimensional interpretation of a 
triangle. As noted in Appendix C, our method dispenses with the need for such 
a heuristic assumption. Similar remarks apply to Barnard's study of curvature. 

Ikeuchi [1983] proposes a method for determining surface orientation for a 
surface that is covered with uniformly repeating texture elements. He assumes 
that the surface element corresponding to tire texture element is planar. Most 
critically it is assumed that the shape of the texture element is known a priori. 
For a particular slant and tilt, the orientation of the projection of the known 
r^ figure changes. Ikeuchi proposes a measure that is superficially similar to the 

symmetry measure in Section 4 to determine tire slant and tilt of a texture 
element. 

Olson [1974] has studied a version of Ames' trapezoidal window illusion. 
Surface slant is judged more accurately when the stimulus is moving consistent 
with the static three-dimensional interpretation than in the purely static case. 
Similarly, Wallach, Weisz, and Adams [1956] have shown that when an ellipse 
is rotated in the image plane about an axis passing through its center and 
normal to the image it is perceived as a spinning oriented circle, like a settling 
spinning penny. The instantaneous surface orientation of the oriented circle 
can be computed by our method. Hildreth's forthcoming PhD thesis on the 
perception of motion discusses the rotating ellipse and proposed a dieoretlcal 
explanation of human perception of it. 
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Appendix A: The maximum likelihood method applied to an ellipse 

Equation (4.30) gives the value F of the logarithm of the Maximum Likelihood 
estimator for surface orientation: 

F = log cos a f dr — f log(cos 2 a -f- sin 2 crcos 2 (a — r))dr. (Al) 

In this Appendix we investigate the slant computed by the Maximum Likelihood 
method for the ellipse 

It is convenient to use the standard parameterization, 
x = a cos 6 

y = bsin6, < < 2ir (A3) 

dr = {a 2 sin 2 + b 2 cos 2 9)*d6 

Pending a more thorough analysis, we assume that the tilt r = tt/2 is computed 
correctly by the Maximum Likelihood method, and restrict attention to slant. 
Witkin [1981, Figure 5] and our own computational experiments suggests that 
this is reasonable; in any case, if it is not so, the inequivalence of the Maximum 
Likelihood and Extremum Principle methods would be assured. Without loss 
of generality, we assume that a > 6, and denote a/b by X > 1. 
Since r = n/2, 



2 1 

cos a = 



1 + tan 2 a 
a 2 sin 2 1 



(AA) 



a 2 sin 2 9 + b 2 cos 2 
since tan a = dy/dx. Also, 



IT 

sin(a — t) = sin(a ) = — cos a, (-A-5) 

and so Eq. (A.l) reduces to 

F = log cos a j dt — j log(l — cos 2 a sin 2 a)dt. ( A.6) 

We are interested in extrema of F, and so we consider dF/da, which we 
write 
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dF £ a 2 sin 2 cos 2 a — b 2 cos 2 , ^ , , A , 

cot a — = f dr, 1 < X (A7) 

da ^ a 2 sin 2 6 cos 2 a + b 2 cos 2 

Suppose that the Maximum Likelihood method is extremized at a = a'. 
The slant that is necessary to interpret the ellipse as a circle is given by 
cos a = 1/X. The Maximum Likelihood method overestimates slant if and 
only if a' > cos —1 (l/X), which is if and only if coser' < 1/X. This is true 
if dFjdc j^ for all a such that 1 < Xcoscr < X. Denote Xcosa by n, 
so that 1 < fj, < X. Substituting in Eq. (A.7) and changing the limits of 
integration, we find 

dF n , r> 2 sin 2 — cos 2 0„ 2 . a „ , i „u ,„ ,. s 

—- = 26/ £- 5 — \ 2 sin 2 + cos 2 <md0 A8 

da ^0 ^2 sin 2 + cos 2 * 

We now split the range of integration into four equal intervals of size tt/A. 
With suitable changes of variables to bring the intervals of integration to [0, f], 
we find 

dF , f% u?cos 2 9 — sin 2 0„ 9 » , . , „,i 

■t- = 4 / ^ : 5- X 2 cos 2 ^ + sin 2 ^ ? 

da Jo M 2 cos 2 + sin 2 , , 

/»"\ /i 2 sin 2 9 — cos 2 „ o . 2 . , 2 .a , . 

' — 5 (X 2 sin 2 0-f-cos 2 0)*d0 

/i 2 sin -j- cos 2 

Over the range of integration, cos 9 > sin 0. Algebraic manipulation of Eq. 
(A.9) shows it to be greater, than zero, so that it is certainly not zero. Hence 
the Maximum Likelihood method overestimates the slant. 



f~\. 



jT\ 



25 



Appendix B: Why square curvature is inappropriate 



In section 3 we rejected the measure 



F 



-j. 



ds 



(5.1) 



because, despite popular belief to the contrary, it does not interpret an ellipse 
as a circle. In this appendix we substantiate that claim. 

Consider the ellipse defined by Eqs. (A.2) and (A.3). The curvature k of the 
ellipse is given by 
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dh 



ds 2 



a 2 b 2 



{a 2 sin 2 6 + b 2 cos 2 0) 3 
where s is the arc length. Substituting Equation (B.3) into (B.l) we have 

a 2 b 2 



(B.2) 
(S.3) 



F(a,b) = J 



2tt 

,'2.™ 2, 



-de 



(BA) 

(o 2 sin z + 6 2 cos 2 0)1 

Suppose, without loss of generality, that the ellipse is in the image plane and 
that a > b. Suppose further that the square curvature performance index (B.l) 
is correct and interprets the ellipse as a circle lying in a plane that is slanted 
with angle a to the y-axis, where 



b = a cos a 



(B.5) 



Equivalently, if we set b' = b/cosa, the measure will choose a such that 
b' = a. Hence if we write b = Xa, and consider F(a, b) to be a function of 
X, the measure should be extremized by X = 1. We will now show that this is 
not the case. 



1 f 2lT 

F(\) = - 

a Jo 



X 2 



■dO 



(sin 2 + X 2 cos 2 0)^ 
Differentiating this with respect to X gives 

r2 " 2X sin 2 9 - 3X 3 cos 2 



(fl.6) 



dF 



-I' 

a Jo 



[sin 2 6 + \ 2 cos 2 6)i 



dO 



(B.l) 
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Evaluating this expression at X = 1 gives 



dF 



x=i 



1 f 2lT 
= - / (2sin 2 0-3cos 2 0)d0 
a Jo 



7T 

a 



(B.8) 



Since the partial of F with respect to X docs not vanish when X is equal to 
one, the circle does not extremize the ellipse for the square curvature measure. 
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Appendix C: The interpretation of some simple shapes 

In this appendix we show that our method correctly interprets a number of 

simple shapes, namely an ellipse, a parallelogram, and a triangle. 

Ellipse 

Suppose the ellipse is given by equations (A.2) and (A.3). It is easy to show 
that the area A and perimeter P arc given by 

A = irab 

P = j 2 \a 2 sin 2 9 + b 2 cos 2 0)4 d* (C,1) 

JO 

Maximizing A/P 2 is equivalent to minimizing Pjs/A. We set b = \a, (X is 
the eccentricity of the ellipse,) and define PjsfA — /(X). We find 

1 /*27r c : n 2 a 

/(X) = — / [-^- + \ cos 2 9)^ d9 (C.2) 

By the same argument presented in appendix B, the ellipse will be interpreted 
as a circle provided X = 1 is a minimum of f(\). Changing die variable of 
the integral to <j> — 9 -}- f we find that 

/(X) = f(l) (C.3) 

which implies that X = 1 extremizes /(X). It is also clear that extrema occur 
in pairs of the form X,l/X. Furthermore, both are stationary points or one 
is a maximum and the other a minimum. Observe that /(X) tends to infinity 
as X tends eidier to zero or infinity. It follows that a sufficient condition for 
X = 1 to be a global minimum is that all pairs of extrema be stationary points. 
Suppose that this is not the case, and let X be the smallest extremum that is 
not a stationary point. Since /(X) tends to infinity as X tends to zero, X must 
be a minimum and 1/Xo a maximum. But 1/Xo is die largest non-stationary 
extremum, and so, by the same argument as above, it must be a minimum. 
This contradiction establishes the result. 
Parallelogram and triangle 

In section 5 we showed that a skew symmetry is always interpreted as an 
oriented symmetry by our method. In particular, a parallelogram is interpreted 
as a rectangle. By die same argument, a rectangle is a skewed symmetry of a 
square. Hence our mediod interprets a parallelogram as a square. 

Similar reasoning shows that a triangle is interpreted as a skewed isosceles 
triangle, which is interpreted as a skewed equilateral triangle. The axes of the 
skewed symmetry join a vertex to the midpoint of the opposite side. Hence 
our method interprets a triangle as an oriented equilateral triangle. 
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Appendix D: The slant and tilt of a skewed symmetry 

In this appendix we calculate the slant a and the tilt r that correspond to the 
oriented real symmetry that is the interpretation of a skewed symmetry whose 
skew angle is 6. 

As a simple, though instructive, example, consider a rhombus of side a and 
included angle 7 (Figure 8). To find the extremizing tilt, we substitute the data 
from Figure 8 into Eq. 4.24, and find 

sin r cost sin(Y — Y)cos(r — 7) 

-I i U i L» = (£».i) 



|l-sin 2 asin 2 r| {l - sin 2 asin 2 (r - 7)}* 



We can rewrite this in the form 



cos 2 r 



<cos 2 (t — 7) — cos 2 r> <cos 2 (r — 7) + 

+ (1 — cos 2 o){\ — cos 2 (t — 7 ))(1 — cos 2 r) — 1 j 

= 

(D.2) 

We assume first that the first factor is zero. It follows that 

cos(r — 7) = i cos t (-D-3) 

and so 

7 — 2r = nir. {DA) 

Since < r < -k and < 7 < tt, there are two possible solutions, namely 

' = \ P.5) 

Observe diat the tilt direction is one of the axes of symmetry o.f the rhombus 
shown in Figure 8. 

Having solved for the tilt, we now solve for the extremizing slant a using 
Eq. 4.23. Recalling from Eq. D.4 that sin 2 (r — 7) = sin 2 r, we find upon 
substitution into Eq. 4.23 

-\\ — sin 2 a sin 2 r> = cos 2 r< 1 — sin 2 a sin 2 r \ (D.7) 
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which we solve to get 

2 cos2 T 



cosher = . (D.8) 



sin 2 r 



The requirement that |cosa| < 1 picks out either Equation (D.5) or (D.6), 
so we get a unique solution. In this case the skew angle is given by 

(D.9) 
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To summarize, 


if 6 < 


we 


get 














T 

coscr 


= 


7T 6 

4 ~ 2 
1 -f sin 5 

cos 5 


and if S > 0, 




















r 


= 


3tt 
4 


<5 
2 








cos a 


— - 


cos 5 



(D.10) 



(D.ll) 

1 + sin 6 ' 

fi These formulae were derived for an equal-sided parallelogram but they will 

clearly apply to the more general case and a rotation through the angles given 
by Equations (D.10) and (D.ll) will unskew any symmetry. It should be noted 
that r is taken to be zero on the axis of symmetry, as in Figure 8. 

To conclude this Appendix, we consider the case that the second factor in 
Eq. (D.2) is zero and die first factor non-zero. We introduce the angle ip by 
analogy with Eq. (2.1) defining orthographic projection: 

tan 2 ib = cos 2 a tan 2 (7 — r) (D.12) 

Then 



. ii x tan 2 ib 

sin (7 — r) = 



cos 2 a -\- tan 2 ib 

o, s cos 2 a 
cos (7 — r) = — — 

cos 2 a -j- tan 1 * ip 



(D.13) 



where, without loss of generality, we suppose < ib < 7r/2. Now by 
assumption, the second factor in Eq. (D.2) is zero: 

cos 2 r -f cos 2 (7 — r) + sin 2 asin 2 rsin 2 (7 — r) — 1 = (-D-14) 
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Figure 8 

A typical skew symmetric figure, namely a rhombus of side a and Included angle 7, with 
skew angle 6 = f — 7. The tilt direction is aligned with one of the axes of symmetry 
of the rhombus, determined by the angle S. 



Using Eq. (D.13), we find 



tan rtan ip = 1, 



from which we deduce 



/"""S 



sin r = cos tp 
cost = fi sin ip, 



(D.15) 



/x = ±l, 

since the ranges of the variables can be assumed to be < a < 7r/2; < 
r < tt;0 < V < tt/2;0 < 7 < tt. From Eq. (D.13) we find that 

. , . , . v tan ib cos a 

sm(7 - r) cos(t — t) = — , f = ±1, 0.16 

cos^ a -f tan" 5 V» 

and, since sin r cos r has the same sign as sin(7 — r) cos(7 — r), //= v. From 
Eq. (D.l) we deduce 



sin T COST 



— sin(7 — r) cos(7 — t) 

1 -sin 2 a sin 2 r* 1 - sin 2 <xsin 2 (7 — t)* 

v sin ^ cos V" 

1 — sin 2 a cos 2 i\r 



(0.17) 



/**% 



f~\ 



/"> 



31 



So far in this Appendix we have only used one of the constraints derived in 
Section 4, namely Eq. (4.24) (from which we derived Eq. (D.l)). We now use 
the second constraint Eq. (4.23), which we can write in the form 



2cos 2 r ,„ . o „ a 
■■ — (1 — sin^a2T)2 

(1 — sin 2 a sin 2 r) i 

2cos 2 (7-r) (£>.18) 



(1 — sin 2 a sin 2 (7 — r)) 2 
— (1 — sin 2 a sin 2 (7 — r))* = 



After some algebraic manipulation, we deduce 



cosa = — tan 2 ip, 



from which it follows that cos a is negative, which is impossible in the range 
under consideration. 
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